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Evidence for the decay X(3872) -^ J/tpuj 
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We present a study of the decays 5"'+ — j> J / ipn^ -k^ n'^ K'^''^ , using 467x10® SB pairs recorded 
with the BABAR detector. We present evidence for the decay mode X(3872) — j- J/V"^, with product 
branching fractions B{B+ -^ X(3872)7^+) x B(X(3872) -> J/ipto) = [0.6 ± 0.2(stat) ± O.l(syst)] x 
10-^ and 6(5" ^ X(3872)is'°) x B(X(3872) ^ J/t/'o;) = [0.6 ± 0.3(stat) ± O.l(syst)] x 10-^ 
A detailed study of the 7r^7r~7r'' mass distribution from X(3872) decay favors a negative-parity 
assignment. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 

The X(3872) meson (denoted in the following as the the possible spin-parity (J^) assignment to 1+ or 2 



X meson) has been observed primarily in its J/ipTT~^TT The decay X — >■ D'^D^w'^ has also been observed 12| 



decay mode [ij-roll • Evidence for its decay to the J/f/'T [ir and interpreted as evidence for X — >■ D*^D'^- this has 

Isl and 'ip{2S)"f [9j final states has established positive C been confirmed by subsequent analyses [ij, llJI- There 

parity. Analyses by the CDF Collaboration of the tt+tt^ has been much theoretical interest in the nature of the 

mass distribution [lOJ , and of the decay angular distribu- X meson |15l422| . Hence, additional experimental infor- 



tion [11] . for the J/^tt+tt decay mode have narrowed mation on new decay modes, especially those sensitive to 



the J^ assignment, is germane to the theoretical under- 
standing of this state. 

In a previous BaBAR publication [2^5, we have con- 
firmed the observation of the y(3940) meson (denoted 
in the following as the Y meson) in the decay mode 
Y — > J /il^Lu reported by the Belle Collaboration in 
-8°'+ -^ Jliiu}K°'+ decay [li]. In the BaBAR analysis 
of this decay mode, the uj — >■ 7r+7r~7r° mass (to3^) region 
was defined as 0.7695 < mg^ < 0.7965 GeV/c^. With this 
requirement and the other selection criteria of Ref. [23| , 
we reported no evidence for the decay X — )■ J/f/'w, al- 
though Monte Carlo (MC) simulation of X-meson de- 
cay to an S'-wave J/tpuj system indicated that this de- 
cay could have been observed. An unpublished Belle 
analysis of B^ — > J / ip-r:^ tt^ tt^ K'^ J7|, which required 
|TO(J/V'7r+7r-7rO) - 3.872| < 0.0165 GeV/c^, reported ev- 
idence for the decay X — )■ J/ijiuj on the basis of 12.4±4.1 
events in the mass interval 0.750 < mg^ < 0.775 GeV/c^. 

In this study we repeat our analysis of the decay modes 



l25j , extending the selected 
tostt region to 0.5 <r7?,37r< 0.9 GeV/c^ in order to investi- 
gate the TO37r distribution in a broader region around the 
Lo meson. 

The data were collected with the BaBAR detector [26| 
at the PEP-II asymmetric-energy e+e^ collider operated 
at the T{AS) resonance. We use the entire integrated 
luminosity at this center-of-mass (cm.) energy (~ 426 
fb~^), which yields a data sample corresponding to about 
467 X 10^ BB pairs. The data were reprocessed with im- 
proved charged-particle-track reconstruction and identi- 
fication. 

The event-selection criteria are identical to those 
in Table I of Ref. [23|, except for the initial ms^ 
requirement. The i3-meson signal region is de- 
fined using the cm. energy difference AS = 
\/s/2, and the beam-energy substituted mass niES 



E*B - 



y/[{s/2+p,-pB)/E,]-^-\pB\ 2 [26|, where (E^p,) is the 



initial state four-momentum vector in the laboratory 
frame (l.f.), \/s is the cm. energy, E'^ is the B meson en- 
ergy in the cm., and ps is its l.f. momentum. Signal B^ 
(B") candidates satisfy |AS| < 20 MeV (15 MeV). In 
events with multiple B candidates (12% of events in the 
region 5.274<mES<5.284 GeV/c^), the candidate with 
the smallest \AE\ is chosen. 

For the i?+-candidate sample, the 7713,1 distribution is 
shown in Fig. [TJ The contribution in each mass inter- 
val is obtained by fitting the corresponding tties dis- 
tribution in the region 5.2 <?7ies< 5.3 GeV/c^ with 
a B~^ signal Gaussian function and an ARGUS back- 
ground function [27[. The Gaussian mean value (/x), 
width {a), and the ARGUS parameter (Carg), are fixed 
to the values obtained when fitting tties for the entire 
J/'07r"'"7r~7r° mass region separately for the 5+ and J5° 
samples (for the B+ sample ^ = 5278.95 ±0.13 MeV/c^, 
a = 2.83 ± 0.14 MeV/c^, and Carg= -37.9 ± 1.8). A 
binned Poisson likelihood fit is performed to the toes 



distribution in each rrig^r interval to obtain the Gaussian 
and ARGUS normalization parameter values, and hence 
to extract the B-meson signal. 

In Fig. [1] there is a small, but clear, ry-meson signal, a 
large w-meson signal, and nothing of significance in be- 
tween. The J /iprj mass distribution shows no significant 
structure, and will not be discussed any further. 
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FIG. 1: The m^^ distribution for B+ -^ J / ^pTT+ n- ■n'^ K+ can- 
didates obtained as described in the text; the dashed vertical 
lines denote the central mass values of the rj- and uj- me- 
son ;28i1. 



In the w-meson region, the signal extends down to 
^ 0.74 GeV/c^; there is also a high- mass tail above 
^ 0.8 GeV/c^, and possibly some small nonresonant con- 
tribution in this region. When we assign oj-Dalitz-plot 
weights to the events in the region 0.74-0.80 GeV/c^, 
the sum of weights (1030±90) is consistent with the signal 
size (1160 ± 60), indicating that any non-w background 
is small, and so we ignore such contributions. Similar 
behavior is observed for 5° decay, but with a selected- 
event sample which is about six times smaller. In the 
following, we define the lower limit of the w-meson mass 
region as 0.74 GeV/c^, but leave the upper limit at 0.7965 
and 0.8055 GeV/c^ for the B+ and B° samples [H, re- 
spectively, in order to focus on the impact of this one 
change on the observed J/V'w mass distribution. The ex- 
tension of the TO3,r region towards lower values increases 
the efficiency slightly. 

The J /ipijj mass distributions for B^-^ — > J/ipcuK'^'^ 
candidates are obtained by using the same fit procedure 
used to obtain the 1713.,^ distribution. We then correct 
the observed signal yields for selection efficiency. Events 
corresponding to B"^+ —^ J/tpujK^'~^ d ecay are created 
by MC simulation, based on Geant4 [30|, in order to 



provide uniform coverage of the the entire rriji^^ range. 
The generated events are subjected to the reconstruction 
and selection procedures apphed to the data. For 5+ 
(5°) decay it is found that the efficiency increases (de- 
creases) gradually from ^ 6% (^ 5%) close to vajj^^ 
threshold to - 7% (- 4%) for mji.^^^ 4.8 GeV/c^. 
Comparison of generated and reconstructed rrij/Tjjuj val- 
ues within each reconstructed ttij/tp^ mass interval en- 
ables the measurement of the rtij j^^ dependence of the 
mass resolution. From a single-Gaussian fit to each dis- 
tribution, the rms deviation is found to degrade gradually 
from 6.5 MeV/c^ at toj/^^-3.84 GeV/c^, to 9 MeV/c^ 
at mj/^„~4.8 GeV/c^. 



The 



iJ/iAi- 



distributions for B'^ 



J/i})ujK^ and 



B^ -^ J/4>ujK^ decay, after efhciency correction in each 
mass interval, are shown in Fig. [UJa) and Fig. [2jb), re- 
spectively. For the latter, corrections for K^ production 
and Kg —> 7r°7r° decay have been incorporated. The 
'mj/ipui range from 3.8425 to 3.9925 GeV/c^ is divided 
into 10 MeV/c^ intervals, while beyond this 50 MeV/c^ 
intervals are used. The same choice of intervals was used 
in Ref. [23|, where the first two were inaccessible, and 
the third was only partly accessible, because of the value 
of the lower limit on TO3^. Clear enhancements are ob- 
served in the vicinity of the X and Y mesons in the 5+ 
distribution, and similar effects are present in the B^ dis- 
tribution, with lower statistical significance. 

The function used to fit the distributions of Fig. [2] 
is a sum of three components. The X meson compo- 
nent is a Gaussian resolution function with fixed rms 
deviation cr = 6.7 MeV/c^ obtained from MC simula- 
tion; the intrinsic width of the X meson (estimated to be 
^ 3 MeV |28j) is ignored. The F-meson intensity con- 
tribution is represented by a relativistic S'-wave Breit- 
Wigner (BW) function [23|. The nonresonant contri- 
bution is described empirically by a Gaussian function 
multiplied by rrij/^p^. The y-meson and nonresonant in- 
tensity contributions are multiplied by the phase space 
factor p X q, where p is the K momentum in the B rest 
frame, and q is the J/^ momentum in the rest frame of 
the J/ip37T system. A simultaneous x^ fit to the distribu- 
tions of Figs. [2ja) and[2jb) is carried out, in which only 
the normalization parameters of the three contributions 
are allowed to differ between Fig. ^a.) and Fig. [2][b). 
The fit describes the data weU (xV^DF = 54.7/51, 
NDF=number of degrees of freedom), as shown by the 
solid curves in Fig. [2] The dashed and dotted curves 
show the X- and K-meson contributions, respectively, 
while the dot-dashed curves represent the nonresonant 
distribution. 

For the X meson, the fitted mass is 3873.0l;i[g(stat) ± 
1.3(syst) MeV/c^, while the mass and width values for 
the Y meson are 3919.lt3:4(stat)±2.0(syst) MeV/c^ and 
3llg°(stat)±5(syst) MeV, respectively. These results are 
consistent with earlier BaBAR measurements p, |23| . 

From the fits of Fig. [21 we obtain product branch- 
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FIG. 2: The corrected rrij/^^ distribution for (a) B^ , (b) i3° 
decays; (c) (inset) shows the low-mass region of (a) in detail. 
The curves indicate the results of the fit. 



ing fraction measurements for 5"'+ -^ XK^'^, X — > 
J I^Lu. The resulting B'^ and B^ product branching frac- 
tion values are [0.6 ± 0.2(stat) ± O.l(syst)] x 10"^, and 
[0.6±0.3(stat)±0.1(syst)] x 10~^ respectively. 

Similarly, we obtain updated values for B{B^ — > 
YK+)xB{Y ^ J/^l)uj) = [3.0l°■^(stat)+o■3(syst)]xl0-^ 
B{B^ -^ YK°) X B(Y -^ J/ituj) = [2.1 ± 0.9(stat) ± 
0.3(syst)] X lO^'^, and for the total {i.e. the sum of 
the X- meson, K-meson, and nonresonant, contributions) 
B{B+ -^ J/iPujK+) = [3.2 ± 0.1(stat)+°-3(syst)] x 10"^ 
and S(B" ^ J/il)uK^) = [2.3 ± 0.3(stat) ± 0.3(syst)Jx 
10^^. These values are consistent with those of Ref. [23J . 
and supersede them. 

Wc define Rx, Ry, and Rnr. as the ratios of the B° to 
5+ branching fractions to the final states XK, YK, and 
nonresonant J/ipujK, and extract these ratios from a si- 
multaneous fit to the data, with the fit function adjusted 
to explicitly contain these parameters. This yields Rx — 



1.0+J];^(stat)+o:2(syst), Ry = 0.7to;;^(stat) ± O.l(syst), 
and Rnr = OJ ± O.l(stat) ± O.l(syst). The values of 
Ry and Rnr are consistent with those in Ref. [23|. The 
statistical uncertainty on Rnr has been reduced signif- 
icantly with respect to Ref. [23| as a result of the in- 
creased luminosity, improvements in event reconstruction 
efficiency, and the use of much larger MC samples in the 
measurement of the selection efficiency as a function of 
mj/,puj, especially for mj/^u> 4 GeV/c^. 

In Ref. @, it was found that B{B+ -^ XK+) x B{X -^ 



J/^JTT'^IT ) 



.5 ± 1.5(stat) ± 0.7(syst)] x 10"^ and 



B{B° ^ XK°) X B{X ^ J/^Pn+TT-) = [3.5± 1.9(stat) ± 



0.4(syst)] X 10~^. We combine these results with those 
from the present analysis to obtain the ratio of the 
branching fractions B{X — > J/ipuj)/B{X -> J/i/)7r+7r~). 
For B+ {B°) events, this ratio is 0.7 ± 0.3 (1.7 ± 1.3), 
where the statistical uncertainties, and those system- 
atic uncertainties which do not cancel in the ratio, 
have been added in quadrature; the weighted average is 
0.8±0.3. This is consistent with that reported in Ref. [7| 
(1.0±0.4(stat)±0.3(syst)). 

In obtaining the quoted systematic errors, systematic 
uncertainties due to tracking (2%), particle identification 
(4.4% and 5.2% for B° and -8+ events), 7r° reconstruction 
efficiency (3.6%), Kg reconstruction efficiency (2%) for 
the B'^ events, and BB event counting (1.1%), have been 
taken into account. The uncertainties on the branch- 
ing fraction values for J/ijj — >■ £~^i~ and a; — >■ 37r 28 1 
have been treated as sources of systematic uncertainty. 
When fitting the thes distributions in each mj/.^^ or mz-„ 
mass interval, the parameters /z, a, and Carg were fixed 
to the values obtained from the fit to the corresponding 
total TTiES distribution. Associated systematic uncertain- 
ties were estimated by increasing and decreasing the cen- 
tral value of each parameter by one standard deviation, 
repeating the analysis, and taking the change in each 
fitted quantity as an estimate of systematic uncertainty. 
Similarly, the systematic uncertainty associated with the 
efficiency-correction procedure was estimated by varying 
its mj/,pi^ dependence within a zblcr envelope, repeating 
the fits to the data of Fig. [21 and taking the corresponding 
changes in fit parameter values as estimates of system- 
atic uncertainty. Additional systematic uncertainties on 
the mass and width of the Y meson were estimated as 
described in Ref. [23J. The main contributions described 
there result from a comparison of the MC input values 
to those obtained after event reconstruction, and from 
the difference in fitted values when a P-wave BW was 
used instead of an S'-wave BW to describe the y-meson 
lineshape. 

Since the AT-moson signal occurs at a low statistical 
level and at very low values of m 



j/^i. 



there is concern 



that the measured signal-event yield might be biased be- 
cause of the low mass tails of the y-meson and nonreso- 
nant contributions. A detailed MC study using samples 
of X-meson events ranging in size from 10-500 events 
showed no evidence of bias, and the spread in extracted 
signal yield was consistent with the corresponding statis- 
tical uncertainty obtained from the fit to the data. 

We now consider the relationship between the X- 
meson signal and the choice of lower mass limit for 
the cj-meson region. In Fig. |3] we show the data cor- 
responding to the first five mass intervals of Fig. [2] 
(3.8425< mj/^^ <3.8925 GeV/c^) before applying the 
efficiency and K^ branching fraction corrections. The 
points shown by open squares indicate the effect of choos- 
ing the 7713^ lower limit to be 0.7695 GeV/c^ rather than 
0.740 GeV/c^. The three lowest intervals then yield no 



signal, and the other two contain only 11 (0.5) events 
in Fig. [HJa) (Fig. ISfb)). This is to be compared with 
42.4]'ly'2 (8.5J^3q) events obtained when the ma^ lower 
limit is 0.74 GeV/c^. Since the number of events in 
Fig. [3] is much smaller than the total number of w-meson 
events (1160 ± 60 for B+ and 206 ± 26 for S° decay), 
and since the rnz-n distribution (Fig. |4ljc)) differs greatly 
from the w-meson lineshape, these might be nonresonant 
37r events. To check the w-meson interpretation, we sum 
the cj-Dalitz-plot weights [29| for the events contribut- 
ing to Fig. [3Ja) (solid points) in the toes signal region 
and obtain 41 ± 13, in good agreement with the number 
from the toes fits. This justifies the w-meson interpreta- 
tion. In contrast, we note that for the 152 ± 20 7y-meson 
events in Fig. [l]the sum of the weights [29| is — 1 ± 42, 
as expected for a uniform Dalitz-plot distribution. 

To determine the significance of the X — j- J /-ipuj signal, 
we extract the signal yields from a fit to the data, prior to 
the corrections for efficiency and K'^ branching fractions, 
as shown in Fig. [3l The fitted values of the masses and 
widths are in agreement with those obtained from the fit 
to the corrected data. An A-meson signal of 21.1 ± 7.0 
events is obtained for B+ decay, and 5.6 ± 3.0 events 
for i3° decay, so that the combined signal is 26.7 ± 7.6 
events. For the combined distribution, the mass region 
3.8625 - 3.8825 GeV/c^ contains 34.0 ± 6.6 events, and 
the fitted curves indicate that only 8.9 ±1.0 events are 
due to the tails of the y-meson and nonresonant distribu- 
tions. We convolve a Gaussian ensemble of background 
Poisson distributions with a Gaussian distribution of ob- 
served events, and obtain probability 3.6xl0~^ that the 
34.0 ±6. 6 events can result from upward background fiuc- 
tuation. This corresponds to a significance of 4.0(7 for a 
normal distribution. On this basis we report evidence for 
the decay mode X — >■ J/^w. 

For the 3.8625 - 3.8825 GeV/c^ region of Fig. [3 we 
plot the TO3,r distributions in Fig. [D Each data point 
results from a fit to the corresponding toes distribution; 
for the points with no error bars, the toes distribution is 
empty. For the combined distribution. Fig. [Jfc), ^ 84% 
of the events have tos^K 0.7695 GeV/c^, the mass hmit 
used in Ref. [23|. The dashed histogram in Fig. SKc) re- 
sults from normalizing the reconstructed A-meson MC 
events to the observed 34 events. Since the J/i/'w sys- 
tem was generated with zero orbital angular momentum, 
this corresponds to positive A-meson parity. One unit 
of orbital angular momentum creates a centrifugal bar- 
rier factor 9^/(1 -I- R^q^) in the description of the J/tpoj 
final state, where i? = 3 GeV~^ is the P-wave Blatt- 
Weisskopf barrier factor radius J31[ (values in the range 
< i? < 5 GeV^^ yield no significant difference). This 
factor suppresses the 7r"'"7r~7r'^ mass spectrum near the 
upper kinematic limit, as shown by the solid histogram of 
Fig- Sic) (also normalized to 34 events). For the dashed 
histogram x^/NDF = 10.17/5 and the x^-distribution 
probability is P(x^, NDF) = 7.1%, while for the solid his- 
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FIG. 3: The uncorrected m 



j/^i. 



mj/^«(GeV/c^) 



distributions for events with 



3.8425 <mj/^^< 3.8925 GeV/c^ for (a) B+ and (b) B° de- 
cays; the open squares correspond to (a) m3,r> 0.7695 (b) 
?Ti3,r> 0.7605 GeV/c'^ [231 . The curves indicate the results of 
the fit. 



togram xVNDF = 3.53/5 and P(x^NDF) = 61.9%. It 
follows that the observed distribution favors the P-wave 
description both quantitatively and qualitatively. If both 
histograms are normalized to the region m^^K 0.7695 



GeV/c^ (which was excluded in Ref. [23|), we expect for 



TO3^> 0.7695 GeV/c^, and hence for the 



^J/4>c 



interval 



3.8725 - 3.8825 GeV/c^, - 4.3 events for the P-wave de- 
scription, and ^ 16.6 events for the S'-wave description. 
However, in Fig. [3] we observe ^ 6 events. In Ref. [32| . 
it was pointed out that for X(3872) -^ D*°D°, the in- 
troduction of one unit of orbital angular momentum in 
the final state could explain the shift in measured X- 
meson mass 1^, |13[ . This observation and the present 
analysis, together with the spin-parity (J^) analysis of 
Ref. fn!], favor J^ = 2" for the X(3872) meson. For 
/ = and J^'^ — 2 ^, the X-meson mass falls within 
the broad range of estimates for the rjrjii^D) charmonium 
state [33J, [SJ] • We conclude that this interpretation is fa- 
vored by the data. 

In summary, we have used the entire BaBAR data sam- 
ple collected at the T(4S') resonance to obtain evidence 
for X — >■ J/tjjuj in B°'^ — >■ J/ipt^K^'^ with product 
branching fraction values [0.6 ± 0.2(stat) ± O.l(syst)] x 
10-5 and [0.6 ± 0.3(stat) ± O.l(syst)] x lO'^ for B+ and 
B^ , respectively. A comparison of the observed mz-n mass 
distribution from X — > J/ifjuj decay to those from MC 
simulations leads us to conclude that the inclusion of one 
unit of orbital angular momentum in the J /tpuj system 
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msTT distribution for events with 
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and (c) the combined distribution. The vertical dashed line 
indicates the nominal a;-meson mass [231 . In (c), the solid 
(dashed) histogram represents reconstructed MC P-wave 
(S'-wave) events normalized to the number of data events. 



significantly improves the description of the data. This in 
turn implies negative par ity for the X meson, and hence 



2 is preferred llj . In addition, we have updated 



the mass and width of the Y meson (3919.ll3'5(stat) ± 
2.0(syst) MeV/c^ and 3ll8°(stat) ± 5(syst) MeV), the 
product branching fraction values for P°'+ — > YK'^-'^ , 
Y — > J/i^uj, and our measurements of the total branch- 
ing fractions for P°'+ — >• J/tpojK'^'^. 
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